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ABSTRACT 

We report observations of Crab giant pulses made with the Australia Telescope Compact Array and 
a baseband recorder system, made simultaneously at two frequencies, 1300 and 1470 MHz. These 
observations were sensitive to pulses with amplitudes 2 3 kJy and widths X 0.5 /is. Our analysis led to 
the detection of more than 700 such bright giant pulses over 3 hours, and using this large sample we 
investigate their amplitude, width, arrival time and energy distributions. The brightest pulse detected in 
our data has a peak amplitude of ~ 45 kJy and a width of ~ 0.5 /is, and therefore an inferred brightness 
temperature of ~ 10 35 K. The duration of giant-pulse emission is typically ~1 /j,s, however it can also be 
as long as 10 fis. The pulse shape at a high time resolution (128 ns) shows rich diversity and complexity 
in structure and is marked by an unusually low degree of scattering. We discuss possible implications 
for scattering due to the nebula, and for underlying structures and electron densities. 

Subject headings: pulsars: general - pulsars: individual (Crab pulsar) - ISM: structure - ISM: 
individual (Crab Nebula) - scattering 



1. INTRODUCTION 

The Crab Nebula pulsar B0531+21 is well known 
for its emission of giant radio pulses, and was origi- 
nally discovered through th e detection of such pulses 
(Staclin fc Reifensteinl 1968h . These sporadic, large- 
amplitude and short-duration bursts can be hundreds or 
thousands of times mo re energetic than regular pulses (e.g. 
iLundgren et al.1 [l995) , and they remain one of the most 
enigmatic aspects of pulsar radio emission. While several 
other pulsars are now known to em i t giant pulses (e.g. 
Coenard et all 119961 : Uohnston etafl 12004 iKnight et all 
20061) . only 3 objects, viz. the Crab and the millisec- 



ond pulsars B1937+21 and J1823-3021A, generate giant 
pulses numerous enough to allow detailed studies of their 
characteristics. 

Observations so far have unravelled several fundamen- 
tal properties of giant pulse emission. It is now fairly 
well established that the fluctuations in their ampli- 
tudes are due to changes in th e coherence of the radio 
emission ()Lundgren et alj Il995l ) , and that they are su- 
perpositions of extre mely narrow nanosecond structures 
(Ha nkins et al.l[2003h . Observations also suggest that the 
emission is broadband, extending over several hundreds 
of MHz (jSallmen et al.lll999t IPopov et alj|2006t ). Another 
important characteristic of giant pulses is their tendency 
to originate in a very narrow phase window of regular radio 
emission; and for the Crab, these windows even coincide 
with the phases of the high-energy (from infrared to 7-ray) 
emission (|Moffett fc Hankinsl [l996). Finally, the distribu- 
tion o f giant pulse energies is known to follow a pow er-law 
form (|Argyle fc Goweij[l972l ILundgren et al.lll995l ). much 
in contrast to the Gaussian or exponential distributio n 
typical of regular pulses (e.g. iHesse fc Wielebinskil[l974h . 

Of all the giant-pulse-emitting pulsars, the Crab is the 
most well-studied. It is known to e mit pulses as ener - 
getic as ~ 10 4 times regular pulses (jCordes et al.l [2004) 
and shows structures persisting down to 2 n s, with inferred 
brightness temperatures as high as 10 37 K (Hank ins et al.l 



l2003h . Given such extreme short durations, giant pulses 
are best studied using data from baseband observations. 
Such data allow coherent dedispersion of voltage samples 
to remove the deleterious effect of interstellar dispersion, 
thereby yielding more accurate descriptions of pulse struc- 
ture, shape and amplitude. However, the bulk of the ob- 
servational studies of the Crab until the early 2000s were 
carried out using traditional filterbank or spectrometer 
data, due to limitations of data throughput and comput- 
ing. With the advent of wide-bandwidth recorders and af- 
fordable high-power computing, these limitations are be- 
ing gradually overcome. As a result, baseband observa- 
tions are increasingly employed in giant pulse s tudies (e.g. 
ISoglasnov et alj|2004t [Popov fc Stappersll2007D . 

In this paper, we report our observations of Crab giant 
pulses made using the Australia Telescope Compact Array 
(ATCA) and the baseband recorder recently developed for 
the Australian Long Baseline Array. We detected over 700 
giant pulses with amplitudes S3 kJy and widths ^0.5 to 
10 /is, and using this large sample we investigate aspects 
such as the pulse amplitude, width, arrival time and en- 
ergy distributions, as well as details of the pulse structure 
and shape. Our observations show a much lower degree 
of scattering than reported before and we discuss possible 
implications for scattering due to the nebula. 

2. OBSERVATIONS AND DATA PROCESSING 

2.1. Observations with the ATCA 

Observations of the Crab pulsar were made in January 
2006. Although the array offers maximal sensitivity for 
pulsar observations in its tied-array mode, our observing 
set-up was designed also to meet certain technical goals 
in addition to pulsar observations, and offered a sensitiv- 
ity equivalent to that of a single ATCA antenna. In this 
set-up, the array was configured to record two indepen- 
dent IF (intermediate frequency) signals, each 32 MHz in 
bandwidth, from a single antenna. At the antenna, two 64- 
MHz bands were Nyquist-sampled at 128 MHz, with four 
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Fig. 1. — Example plots illustrating our giant pulse detection. 
Top: Dedispersed time series of a short data segment around the 
giant pulse. Middle panels: Diagnostic plots from the giant pulse 
search; plots of signal-to-noise ratio vs trial dispersion measure and 
time resolution, and dedispersed pulses over four 8-MHz sub-bands 
of the 32 MHz recorder bandwidth. Bottom: The dedispersed pulse 
over a time window of 33 ms (i.e. one rotation period). 

bits per voltage sample, for two orthogonal linear polarisa- 
tions. These data were converted to circular polarisations 
by inserting a 90° phase shifter in the signal path. The sig- 
nals were then digitally filtered to provide two-bit sampled 
voltage time series for the 32 MHz bandwidth. This results 
in an aggregate data rate of 512 megabits per second. The 
system provided flexibility to record either two separate 
frequency bands from a single antenna, or two identical 
frequency bands from two different antennas. The pulsar 
observations were made over a total duration of 3 hours; 
the first half of which used two 32-MHz bands (dual polar- 
isation) centred at frequencies 1300 and 1470 MHz, while 
the remaining half used two identical frequency bands cen- 
tred at 1300 MHz from two different antennas. This ac- 
cumulated a total data volume of 5.5 terabytes, but split 
into a large number of short blocks of 10 s duration each to 
facilitate further analysis and processing. The data were 
stored onto disks for offline processing and all processing 



operations, such as coherent dedispersion, detection and 
search for giant pulses, were performed at the Swinburne 
University of Technology Supercomputing Facility. 

2.2. Data Decoding and Dedispersion 

Data samples from both the polarisation channels of the 
two IF bands were organised into a single data stream 
while recording, and were subsequently decoded into sep- 
arate IF data streams before de-dispersion and further pro- 
cessing. As our processing software is designed to handle 
a single data stream at a given time, this meant two sep- 
arate passes for processing each short block of data. The 
unpacked data from this stage were then dedispersed using 
the cohe rent dispersion removal technique originally devel- 
oped by Ha nkinsl (|l97lh . The pr ocedures we have ado pted 
are similar to tho se described in iKnight et alj (|2006l ) and 
iBhat et "all (J20O7). Voltage samples are first Fourier trans- 
formed to the frequency domain and the spectra are di- 
vided into a series of sub-bands. Each sub-band is then 
mult iplied by an inverse response filter (kernel) f or the ISM 
(e.g. lHankins k Rickettlll975Hvan Stratenll2003h . and then 
Fourier transformed back to the time domain to construct 
a time series with a time resolution coarser than the orig- 
inal data. By splitting the input signal into several sub- 
bands (four in our case), the dispersive smearing is essen- 
tially reduced to that of an individual sub-band. This also 
means the procedure uses shorter transforms than single- 
channel coherent dedispersion and is consequently a more 
computationally efficient method. The coherent filterbank 
stream data obtained in this manner are then square-law 
detected, corrected for dispersive delays between the sub- 
bands, and summed in polarisation to construct a single 
coherently dedispersed time series for the entire band. 

2.3. Search for Giant Pulses 

Following the procedures of decoding and dedispersion, 
we performed a rigorous search for giant pulses within each 
10 s block of data. Our pulse detection procedure involved 
progressive smoothing of time series with matched filters 
of widths ranging from 0.5 to 16 /xs in steps of 0.5 /xs, 
and identifying the intensity samples that exceed a set 
threshold (e.g. 12<r for the 0.5 /is smoothing time). In 
addition to performing dedispersion at the Crab's nominal 
DM of 56.75 pc cm -3 , we also performed this procedure 
over a large number of adjacent DM values (typically over 
a DM range ~ 0.5 pc cm~ 3 , in steps of 0.001 pc cm -3 ). 
For each DM and the matched filter width, we computed 
the signal-to- noise ratio (S/N) of the pulse amplitude over 
a short stretch of data centred at the pulse maximum. 
From this analysis, diagnostic plots are generated for each 
candidate giant pulse as shown in Fig. [T] These plots are 
subjected to a careful human scrutiny to discriminate real 
giant pulses from spurious signals. 

2.4. Summary of Detections 

We adopted a rather conservative threshold of 12cr in 
order to account for possible departures of noise statis- 
tics from those expected of pure white noise and also to 
limit the number of candidate signals to a reasonable num- 
ber. With this threshold and the 0.5 fis final time reso- 
lution (At) employed in our analysis, a pulse will need 
to be S3 kJy in amplitude to enable a clear detection 
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(see § 12. 5p . We detected 706 giant pulses from our 3-hr 
long observation, of which 413 are from the first half of 
observations that were made simultaneously at 1300 and 
1470 MHz. However, only 70% of these 413 were detected 
in both freq uency bands. This fr action is similar to that 
reported bv lSallmen et al.l (|l999f ) from their observations 
at 600 and 1400 MHz, although quite different fractions 
have been r eported from observa tions at widely separated 
frequencies (iKostyuk et al. l l2003h and at lower frequencies 
(|Popov et al.ll2006h . This probably suggests that not all 
pulses are entirely broadband, a characteristic presumably 
intrinsic to the giant-pulse emission phenomenon. Our 
sample includes pulses with widths ranging from 0.5 to 10 
[is. Given At = 0.5 [is, any shorter-duration pulses that 
may be present in our data will be smoothed to an effective 
resolution^ of 0.5 or 1 /is. Moreover, pulses broader than 
10 us are likely to be weaker than our detection thresh- 
old. The strongest pulse in our observations is 0.5 /xs in 
duration and has a S/N of 217. 

2.5. System Sensitivity and Flux Calibration 

The Crab Nebula is a fairly bright and extended source 
in the radio sky, with a flux density of ^955 v~ - 27 Jy 
(Bictcnh olz et al.l (|l997f ); where v is the frequency in GHz) 
and a characteristic diameter of ~5'.5. Thus, in general, 
there will be a significant contribution from the nebula to 
the system noise (S), depending on the frequency of obser- 
vation and the coverage of the nebula within the telescope 
beam. However, in our case, the problem is much simpli- 
fied as the nebula is unresolved by a single antenna of the 
array (half-power beam width « 33' at our observing fre- 
quencies). In fact, the entire nebular region occupies only 
a small fraction (« 3%) of the beam solid angle. Measure- 
ments made in parallel with the observations yield system 
temperatures (T sys ) of 114 ± 3 and 103 ± 2 K respec- 
tively at 1300 and 1470 MHz. Thus, assuming a nominal 
gain (G) of « 0.1 K Jy^ 1 for a single ATCA antenna at L- 
band, these measurements translate to system equivalent 
flux densities (S sys =T sys /G) of 1140 and 1030 Jy respec- 
tively at 1300 and 1470 MHz. Scaling for our processing 
parameters and the recording bandwidth (AB), and also 
accounting for the loss of S/N due to our 2-bit digitisation, 
these estimates will correspond to system noise of 250 and 
227 Jy at the l-er levefl In other words, minimum de- 
tectable pulse amplitudes of 3 and 2.7 kJy respectively at 
1300 and 1470 MHz for a 12-cr threshold. 

3. STATISTICAL PROPERTIES OF GIANT PULSES 

3.1. Pulse Amplitudes and Energies 

An important distinguishing characteristic of giant 
pulses (GPs) is their power-law distribution of ampli- 
tudes and energies. The amplitude determination can be 
potentially influenced by factors such as smearing due to 
instrument and any residual dispersion, as well as external 
effects such as scattering and scintillation due to the ISM. 
i 

In practice, pulses that are narrower than 0.5 fis will be detected 
as either 0.5 fis wide (most flux in one sample) or 1 /is wide (pulse 
comes in just on the border of two samples) in our analysis. 
2 

S = rjN a T aya G _1 (n po! ABAt) -1 / 2 , where r? is the loss of S/N due 
to digitisation, is the detection threshold in units of a, n po i is 
the number of polarisations summed. 
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Fig. 2. — Histograms of giant-pulse amplitudes at 1300 and 1470 
MHz (in red and green respectively) constructed from our sample 
of giant pulses. Estimates for the best fit slope (/3) are —2.33 ±0.14 
and -2.20±0.18 at 1300 and 1470 MHz respectively, and are shown 
as short-dashed and long-dashed lines. The brightest pulse detected 
in our data (at 1300 MHz) has a peak amplitude of ~45 kJy. 

As a result, oftentimes measurements te nd to underesti- 
mate t he t rue amplitudes. Obs ervations of lLundgren et al.l 
(|!995h and lCordes et al.l (|2004l ) (at 430 MHz) were marked 
by significant instrumental and dispersion smearing (315 
and 152 [is respectively) and were more prone to scintil- 
lation and scattering given their observing frequencies S 1 
GHz. On the other hand, smearing due to instrument or 
residual dispersion is negligibly small for baseband obser- 
vations, and consequently our data enable a better and 
more accurate characterisation of the pulse amplitudes. 
Fig. [2] shows the distributions of pulse amplitudes for our 
data. 

The cumulative energy distribution provides a mean- 
ingful way of ch aracterizing the freq uency of occurrence of 
GPs. Following iKnight et all ((2006) . we define this distri- 
bution in terms of the probability of a pulse having energy 
greater than Eq , and can be expressed as 

P(E>E Q ) = KE% (1) 

Figure [3] shows such a distribution for our data. The pulse 
energy is estimated by integrating the amplitude bins over 
the extent of emission. The uncertainties in the energy 
estimates depend on the pulse strength and width, and 
may range from ~0.5% for strong and narrow pulses, to 
as much as ^35% for weak, broad pulses. Such large uncer- 
tainties at low energies, and possible modulations in pulse 
amplitudes due to scintillation, may probably explain a 
gradual flattening seen at low energies. In any case, much 
in agreement with the earlier work, no evidence is seen for 
a high-energy cut-off, suggesting that exceedingly bright 
and energetic pulses are potentially observable over longer 
durations of observation. 

Excluding the segments near the low and high ends of 
the distribution where the behaviour tends to depart from 
a power-law, we obtain a best-fit value of —1.59 ± 0.01 
for a (where the error is purely formal) and K — 3.97. 
Interestingly, this slope is comparable to the published 
values for other GP-emitters such as PSRs J0218+4232 
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Fig. 3. — Cumulative distribution of giant-pulse energies at 1300 
MHz. The slope estimates of the linear segments below and above 
an apparent break near ~10 kjy fis are —1.4 and —1.9 respectively, 
while the overall slope is —1.6. 



(iKnight et al.ll2006T). B0540-69 (I Johnston et alJl2004ft a nd 
B1937+21 (jCognard et aHll996t ISoglasnov et al.ll2004ft 

However, it is known from the work of iMoffettJ (| 19971 ) 
that a single slope does not accurately describe the energy 
distribution at 1400 MHz and that there occurs a break 
around an energy of ^2 kjy fis, where the slope a changes 
from —3 to —1.8. Given that this energy corresponds to 
the detection threshold of our GP search (see § I2.5[) , our 
slope estimate of —1.6 is i n general agreement with Mof- 
fett's work. Recent work of lPopov fc Stappersl (|2007f ) sug- 
gests that the slope estimate depends on the pulse width, 
evolving from —1.7 to —3.2 when going from their shortest 
(4 fis) to longest (65 fis) GPs. Their analysis also shows 
a clear evidence for a break where the power-law tends 
to steepen (their Fig. 2). Such a break is also apparent 
in our data, as the slope tends to steepen near ~10 kjy 
fis. Adopting this as the break point, we obtain slope esti- 
mates of -1.37 ± 0.01 and -1.88 ± 0.02 for the segments 
below and above it (K = 3.72 and 4.42 respectively). 

The pulse energy distribution also serves as a useful 
guide for estimating the rates of occurrence of giant pulses. 
For instance, on the basis of our observations, a pulse that 
is about 50 times more energetic than regular pulsed can 
bo expected roughly once in 8500 pulse periods, i.e., ap- 
proximately once every 5 minutes. The most energetic 
pulse detected in our data has an estimated energy of 66 
kjy fis; i.e. an energy per fis that is almost ~ 10 5 times 
larger than that of regular pulses. 

The brightest pulse detected in our data has an esti- 
mated peak flux density (S v ) of ~ 45 kjy and an effective 
width of 0.5 (is. The equivalent brightness temperature 
(Tb ) is given by (based on the light-travel size and ignor- 
ing relativistic dilation), 



2k T 



where v is the frequency of observation, St is the pulse 
width, ke is the Boltzmann constant and D is the Earth- 
pulsar distance. The inferred brightness temperature for 
our strongest pulse is therefore ~ 10 35 K. This is almost 

3 

Assuming mean pu lse energies estimated from parameters in the 
published literature (Manchester et al. 20Q5|). 



~1000 times brighter than the brightest pu l se det ected 
in the Arecibo observations of ICordes et al.l (|2004h . To 
the best of our knowledge, this marks the brightest pulse 
ever recorded from the Crab pulsar within the L-band fre- 
quency range (i.e. ~l-2 GHz). Detection of such exces- 
sively bright giant pulses offers a pr omising technique for 
finding pulsars in extern al galaxies (|Johnston fc Romanil 
2003; ICordes et al.ll2004ft . 



3.2. Pulse Widths 

There have been very few attempts o f characterising, 
iant-pulse widths. The first analysis bv iLundgren et all 
1995ft was limited by their coarse time sampling (205 
fis) and was further constrained by a severe dispersion 
smear ing (70 fis). More recent work of iPopov fc Stappersl 
(|2007ft had a much higher time resolution (4.1 fis). As 
our observations were made at a higher frequency than 
the above (hence the effect of scattering is reduced) and a 
final time resolution (0.5 fis) that is eight-fold better than 
IPopov fc Stappersl (|2007f ). our analysis allows a more ac- 
curate characterization of giant-pulse widths. 

Our observations show that many giant pulses tend to 
have a significant structure at 0.5 fis resolution, ranging 
from a simple narrow spike to seve ral closely spaced com- 
pon ents within a few fis. Follo wing [Lundgren et al.l |l995) 
and lPopov fc Stappersl (|2007l ). we adopt the notion of "ef- 
fective pulse width" W e , which is essentially the averaging 
time that yields the maximum S/N. In most cases, it is a 
close representation of true pulse width, although it may 
be slightly biased to stronger component (s) in the case 
of highly structured pulses. The histogram of widths ob- 
tained in this manner is shown in Fig. [¥] The distribution 
is highly skewed, with measured widths ranging from 0.5 
to 10 fis and with a clear peak at 1 fis. Giant pulses of 
widths larger than 10 fis are absent in our data; however, 
this may be a selection effect as such pulses may very well 
be below our detection threshold. Our analysis confirms 
a general tendency for str onger pulses to be narrow er, an 
observation also noted by lPopov fc Stappersl (|2007l ). 

Our distribution of giant-pulse widths can be compared 
with a similar plot obtaine d for GPs from the mi llisecond 
pulsar PSR B1937+21 bv ISoglasnov et all (j2004D. Simi- 
larities in the two distributions are quite striking, despite 
the fact that GPs from PSR B1937+21 are intrinsically 
narrower (W e ~ 15 ns) and show little structure. Thus, it 
appears that exponential-tailed distributions are probably 




Fig. 4. — Left: Distributions of giant-pulse widths at 1300 and 
1470 MHz (in red and green respectively). Right: Plot of pulse 
amplitudes vs the effective pulse widths. An apparent low cut-off 
at 0.5 fis and a discretization in pulse width are due to the 0.5 /is 
final time resolution adopted for our processing. 
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common characteristics of giant-pulse widths. Studies of 
more pulsars are necessary to confirm such a conjecture. 

3.3. Arrival Times and Phases 

One of the most striking properties of Crab GPs is their 
occurrence within well-defined, relatively narrow longitude 
ranges of regular radio emission. This is also true with the 
millisecond pulsars PSRs B1937+21 and J0218+4232, ex- 
cept that in the case of the former, GPs tend to occur 
near the trailing edges of the main pulse and interpulsc 
(Sog lasnov et al.|[2004D . and in the latter case, they oc- 
cur near the two m inima of the integrated pulse profile 
(Knight ct al. 200^). In order to determine the pulse ar- 
rival times, we use the pulsar's spin-down model along with 
the TEMPO software package to obtain the anticipated pulse 
arrival phases. As most pulses in our data are very nar- 
row, our data allow precise determinations of the arrival 
times and phases. Pulse times-of-arrival (TOAs) obtained 
in this manner are plotted in Fig. 

The total range of longitudes where GPs occur (in the 
main pulse region) is approximately ± 200 /is, or ±2°. 2 in 
angular units, with an RMS of 84 /xs (0°.9). Thus, a vast 
majority of GPs (75%) occur within a narrower window 
of ±100 fj,s. As contributions from the dispersive and in- 
strumental smearing are negligible in our case, this RMS 
corresponds essentially to the intrinsic pulse-phase jitter 
and can be compared to a value of 90 /is estimated by 
iLundgren et all (|l995l ) for their data at 800 MHz. The plot 
of joint statistics of timing residuals and pulse amplitudes 
shows no obvious correlation except for a general tendency 
for stronger pulses to originate within narro wer phase win- 
dows. A similar property was also noted by ICordes et al.l 
( 2004]) in their Arecibo observations at a frequency of 430 
MHz. Finally, a vast majority of GPs tend to occur within 
the main pulse window — 87% at the main pulse region 
and the remainder 13% in the interpulse region. This is in 
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Fig. 5. — Arrival times of giant pulses detected in the main pulse 
and interpulse regions (the red and green symbols respectively) at 
1300 MHz (left), along with their statistics (right, top) and joint 
statistics with pulse amplitudes (right, bottom). 
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Fig. 6. — A bright giant pulse at 1300 and 1470 MHz, where the 
plotted time window corresponds to RiO.003 cycles of the pulsar's 
rotation period. This pulse pair is displayed at a time resolution of 
128 ns and illustrates significant differences seen in pulse width and 
structural details between the two frequencies. 



excellent agreement with the pub lished values in the lit- 
erature (|Popov fc Stapperd (120071): 84% a nd 16% respec- 
tively at 1200 MHz; iKostvuk et~atl (|2003h : 86% and 14% 
respectively at 2228 MHz). 



4. PULSE BROADENING AND SCATTERING DUE TO THE 
NEBULAR PLASMA 

4.1. Pulse Shape and Estimation of Scattering 

Figure [6] shows an example of a bright giant pulse from 
our data. At a time resolution of 128 ns, the pulse is 
resolved into multiple narrow components, and this readily 
confirms the basic picture of GPs comprising fine structure 
on very short timescales. Discerning such fine structure is 
however limited by time resolution and the smearing due 
to multipath scattering. This pulse pair also exemplifies 
the differences seen in pulse structure between 1300 and 
1470 MHz, a detailed analysis of which is beyond the scope 
of the present wo rk. Our pulse shapes can be compared 
to observations of ISallmen et al.l |l999) ~10 yr ago (their 
Fig. 1), and it is striking that our observations are marked 
by a much lower degree of scattering. 

In order to estimate the pulse broadening due to inter- 
stella r scattering (e.g. IWilliamsonl Il972t ICordes fc Laziol 
l200lh . we adopt th e CLEAN-based deconvolution ap- 
proach developed bv lBhat et al.l (|2003l ). Unl ike the tradi- 
tional frequency-extrapol ation approach (e.g. lLohmer et all 
[200l [K uzmin et al.ll2002t ), this method makes no prior as- 
sumption of the intrinsic pulse shape, and thus offers a 
more robust means of determining the underlying pulse 
broadening function (PBF). The procedure involves de- 
convolving the measured pulse shape in a manner quite 
similar to the CLEAN algorithm used in synthesis imag- 
ing, while searching for the best-fit PBF and recovering 
the intrinsic pulse shape. It relics on a set of figures of 
merit that are defined in terms of positivity and symmetry 
of the resultant deconvolved pulse and some parameters 
characterizing the noise statistics in order to determine 
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Fig. 7. — A bright giant pulse from observations at 1300 MHz. 
The time window corresponds to 0.002 cycles of the pulse period. 
The reconstructed pulse (green) and the best-fit pulse-broadening 
function (inset) from our deconvolution procedure are also shown. 
The e _1 point of this one-sided exponential gives an estimate of the 
pulse-broadening time, for which our measured value is 0.8 ±0.4 /is. 



the best- fit PBF. For the purpose of our analysis, we as- 
sume the simplest and most commonly used form for the 
PBF that corresponds to a thin-slab scattering screen ge- 
ometry The func tional form for su ch a PBF is a one-sided 
exponential (e.g. I Williamson! 1 1 9721 ) and is given by 



au) = [ ±-) ex P n/i. 



(3) 



where U(t) is the unit step function, U(< 0) = 0, U(> 
0) = 1. Fig. [7] shows the best-fit PBF and the recon- 
structed pulse obtained in this manner. The e _1 point of 
this PBF is the scattering time , for which our decon- 
volution procedure yields an estimate of 0.8 ± 0.4 fxs. 

This measurement of a low scattering is further sup- 
ported by the following observational facts. First, a pulse 
broadening of ~0.8 /is will imply rapid intensity decorrela- 
tions in frequency on characteristic scales ~ l/27rrd ~ 0.1 
to 0.3 MHz and this is confirmed by our analysis. Second, 
as seen in Fig. [H a large number of pulses detected in our 
data have widths ~1 fjs, just about what we would expect 
given the ~0.8 /is broadening due to scattering and the 
0.5 fj,s time resolution. 

4.2. Scattering due to the Nebula 

Our me asurement of th can be compared to that re- 
ported by ISallmen et alJ (| 19991 ) based on their observa- 
tions in 1996. Their value of 95 ± 5 ^s at 600 MHz 
would scale to 4.3 fj,s at 1300 MHz, assuming a canon- 
ical dependence expected from scattering due to a 
turbulent plasma screen. This is five times larger than our 
measurement. A somewhat smaller value i s expe cted on 
the basis of observations of iKuzmin et alJ (|2002l ) in the 
year 20 00. However, our measurement is consistent with 
that of iBhat et all (|2007l) at 200 MHz, if we extrapolate 
their measurements using their revised frequency scaling 

Of z/-3 5±0.2 

As described in § 12.31 our GP detection procedure in- 
cludes determination of the best DM by performing dedis- 
persion over many trial values around its nominal value. 
As most pulses in our data are very narrow, this proce- 
dure allows a precise determination of the Crab's true DM 
at our observing epoch. An estimate of 56.751 ± 0.001 
pc cm -3 obtained in this manner is further confirmed by 



measuring the time delay between the pulse arrival times 
at 1300 and 1470 MHz. This value is in excellent agree- 
ment with that reporte d in the Jodrell Bank (JB) monthly 
ephemeris (|Lvnd fl982) on the nearest date of our observ- 
ing, but it is significant ly lower than that m easured near 
the observing epochs of ISallmen et al.l (jl999l ). 

Thus, over a 10-yr time span between Sallmcn et al's 
and our observations, the pulsar's DM decreased by 0.09 
pc cm" 3 while the pulse broadening decreased by a fac- 
tor of 5. Given the strong frequency dependence of pulse 
broadening (t<j oc v~ x ) and the observation al evidence for 
the s c aling index :r cha nging over the time (jKuzmin et al.l 
I2002t IBhat et all 12007ft . it is more meaningful to adopt a 
frequency independent parameter such as the scattering 
measure (SM) for comparison purposes. This parameter 
quantifies the total scattering along the line of sight (LOS) 
and is defined as the LOS integral of C 2 , which is the 
spectral coefficient of the wavenumber spectrum of elec- 
tron density irregularities. It can be related to Td via the 
relation r d » 1.1 W T SM 6/5 v- x D, where v is in GHz, D is 
in kpc, and W T is a geometric factor th at depends on the 
LOS- distribution of scattering material (jCordes fc Rickettl 
Il998ft . Assuming W T — 1, we can estimate the effective 
SM for a uniform medium, and for the measurements of 
Sallmen et al's and ours, we obtain values of 1.3 x 10~ 2 
and 3 x 10~ 3 kpc m~ 20 / 3 respectively. That is, the SM 
changed by 0.01 kpc m -20 / 3 when the DM changed by 0.09 
pc cm" 3 , or a factor 4 decrease in the scattering strength 
associated with a DM change of only 0.16%. 

A closer look at the JB ephemeris reveals quite a system- 
atic decrease in the Crab's DM variations between 1996 
and 2006, with the lowest DM recorded near mid-2005 and 
a reversal of the trend in early 2006. A plot of SM esti- 
mates for the available measurements during this period 
is shown in Fig. [8] along with the DM measurements at rel- 
evant epochs. These plots illustrate a gradual reduction in 
the Crab's DM and scattering over this 10-yr time span. 
Such variations are too large and smooth to be caused by 
refractive scintillation effects in the ISM, and are therefore 
indicative of a nebular origin. 

The material within the Crab Nebula, specifically the 
perturbed thermal plasma associated with it, has often 
been advanced as the source of excessive scattering and 
anomalous DM variations on several occasions. The most 
remarkable observations in support of thi s are the anoma- 
lous scattering recorded in 19 74-1975 (|Lvne fc Thornd 
Il975t llsaacman fc Rankin 19771) and the refl ection event 
in 1997 (|Backer et al.ll2000t iLvne et ai1l200lft . The scat- 
tering event in 1974 was especially noted for its extreme 
activity, where the pulsar's DM rose by 0.07 pc cm" 3 and 
the scattering increased by an order an magnitude over a 
time span of several months, llsaacman fc Rankinl (j!977ft 
ascribed this to a two-component scattering-screen model 
where the variable screen associated with the nebula gives 
rise to such rapid changes in both DM and scattering. The 
anomalous dispersion event in 1997 was seen as discrete 
moving echoes of the pulse, and was interpreted as reflec- 
tio ns from an ionized shell in the outer parts of the nebula 
by ILvne et "all (|200lh , and in terms of variable optics of a 
triangular prism located in the interface of t he nebula and 
the supernova ejecta by I Backer et alJ (|2000t ). 

Fig. [S] reveals a similar change in DM and a large but 
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Fig. 8. — Measurements of DM and SM of the Crab pul- 
sar from observations over the past 10 years (1996 to 2006). 
The SM values (upper panel) are estimated from the published 
pulse-br oadening measurem ents and dccorrclation bandw idths - 
1996.4: [Sa llmcn ct afl fl998h : 2000.6: IKostvuk ct al. (200.3); 2001.5: 
I Kuzmin et all J200^) ; 2002.1 2002.4: ICordes ct al. (2004); 2003.9: 
IPqpov et al.l 120061) ; 2005.6: IBhat et al.l (120071) : 2006.0: our obser- 
vations. The DM measurements are from the Jodrell Bank monthly 
ephemeris (lower panel); the thick solid circles correspond to the 
observing epochs of the SM measurements. 



less dramatic change in scattering. A direct interpretation 
of such observations can be made along the lines of a vari- 
able scattering screen model, as we discuss in detail in the 
following section. 

4.3. Implications for the Nebular Structure and Densities 

The observed changes in DM and SM can be used to 
constrain the combination of the electron density and the 
"fluctuation parameter" in the nebular region, denoted 
as n e cn and F cn respec t ively. Following ICordes fe Laziol 
(|2003h and IBhat et alj (|2004D , the measured decrements 
in DM and SM can be expressed as SDM — n e ^ cn 8s 
and SSM = C% cn 5s , where Ss is the size of the neb- 
ular scattering region and C 2 xn is the equivalent turbu- 
lent intensity (assuming a uniform distribution of scat- 



tering mat erial). The parameter s F cn and C 2 cn can be 
related as (|Tavlor fc Corded 119931 : ICordes fc Laziol l2002t ) 
Cl, cn = C SM F cn nl cn , where C SM = [3(2^) 1 /3]-i ^ 
for a Kolmogorov spectrum, and K u =10.2 m _20 / 3 cm 6 to 
yield SM in units of kpc m -20 / 3 . The above expressions 
can be combined to yield the ratio of SSM and SDM , and 
is given by 

Thus our measurements of 5SM= 0.01 kpc m~ 20 / 3 and 
SDM = 0.09 pc cm -3 yield F cn n ejCn = 55. The fluctuation 



parameter is essentially a product of normalized variances 
(at small and large scales) and other terms such as the 
outer scale and filling factor. The electron density n e c „ is 
unknown, but assuming a nominal value of 1 cm -3 , we 
get F cn = 55, which is much larger than that typical of the 
Galactic spiral arms (F ~ 10). 

The large values of SSM and F cn indeed confirm the ma- 
terial within the nebula as the source of excessively strong 
scattering. While the scattering and dispersion events of 
1974 and 1997 were interpreted in terms of a single large 
structure with electron density ^1,500 cm -3 , the long- 
term systematic variations in dispersion and scattering as 
shown in Fig. [5] can be interpreted in terms of a scenario 
whereby the nebular segment of the LOS is populated by 
many smaller structures of much lower densities. A varia- 
tion in the number density of such structures may then ac- 
count for the observed changes in DM and SM. As the neb- 
ula is thought to comprise fine structure on many length 
scales, perhaps even on scales much finer than the fila- 
mentary structure suggested by optical observations (e.g. 
iHester et al.lll995h . such a picture seems quite plausible. 

For the sake of simplicity, if we model the measured 
changes in scattering and dispersion to arise from N such 
structures of size i5s sc and density n ejSC , the resultant con- 
tributions to SM and DM are given by 



=N 



5SM=y SM SC i — N Csm F sc n 



i=N 



SDM = 



DM sr ,.i =Nn e 



, Ss„ 



(5) 



(6) 



Following equations (4)-(6) and the constraint F cn n e ^ cn = 
55, the electron density of such a structure can be esti- 
mated as 

( SSM \( 1 \ (?) 



c ^ V55 C SM J \N5s s 
Thus, with just two direct measurements alone (SDM and 
SSM ) , it is hard to constrain all 3 free parameters of the 
model. However, assuming a reasonable value for Ss sc ~ 
10 -5 pc (i.e. an order of magnitude smaller than that 
implied by the 1997 reflection event), we estimate n ejSC ~ 
100 cm~ 3 for N^IOO. This is almost an order of magnitude 
smaller than the densities required to produce the reflec- 
tion event. Indeed, several different combinations of size 
and density are possible; nonetheless, the underlying pic- 
ture is the presence of many moderately-dense structures 
in the nebular region, with a filling factor ~ 10 -3 . 

We note that the electron density estimated above, es- 
timates of the elect r on temperature in the Crab nebula 
(|Temim et al.l 120061 : IHester et al.l 1 19951 iDavidsonl fl97l 
that yield values between 6000 and 16000 K, and the size of 
the Crab nebula (approximately 6' = 3 pc, for a distance of 
2 kpc), imply that the nebula should produce a minimal 
optical depth to free-free absorption at GHz frequencies 
and a significant optical depth at lower frequencies. We 
estimate, based on these parameters, that a free-free opti- 
cal depth (t//) of 0.007 at 1 GHz may be possible, along 
lines of sight that include the scattering structures dis- 
cussed above. For radio sources behind the Crab nebula, 
this would give a decrease in the observed flux density, 
compared to the intrinsic flux density, of 0.7% due to free- 
free absorption at this frequency. At lower frequencies this 
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decrease would be more significant: 3% at 500 MHz, 12% 
at 250 MHz, and 60% at 100 MHz. 

If a radio interferometer operating at these frequencies 
could achieve an angular resolution high enough to resolve 
the Crab nebula and detect compact sources of emission 
behind the nebula, it may be possible to survey the free- 
free absorption due to the nebula along many lines of sight, 
probing the structures that are producing the scattering 
of the pulsar emission. The upcoming future instruments 
such as the extended LOFAR telescope in The Netherlands 
and Europe, or the Murchison Wide-field Array (MWA) in 
Western Australia may be able to undertake such a survey. 

5. SUMMARY AND CONCLUSIONS 

Using the ATCA and a baseband recording system, we 
detected more than 700 giant pulses from the Crab pul- 
sar from our continuous and uniform recording at 1300 and 
1470 MHz over 3 hours. This large sample is used for inves- 
tigating statistical properties of giant pulses, such as their 
amplitude, width, arrival time and energy distributions. 
The amplitude distribution follows roughly a power-law 
with a slope of —2.33 ± 0.15, which is shallower compared 
to those from previous observations at lower frequencies. 
The pulse widths show an exponential-tailed distribution 
and there is a tendency for stronger pulses to be narrower. 
A majority of pulses (87%) tend to occur within a narrow 
phase window (±200 fis) of the main-pulse region. Finally, 
the distribution of pulse energies follows a power-law with 
a slope of —1.6, and there is evidence for a break near ~10 
kJy (is. 

The brightest pulse detected in our data has a peak am- 
plitude of 45 kJy and a width of 0.5 /is, implying a bright- 
ness temperature of 10 35 K, which makes it the brightest 
pulse recorded from the Crab pulsar at the L-band fre- 
quencies (1-2 GHz). 

Our observations show that many giant pulses comprise 
multiple narrow components at a time resolution ~128 ns, 
which confirms the fundamental picture of giant pulses 
being superpositions of extremely narrow bursts. Fur- 
ther, the measured pulse shape is marked by an unusually 
low degree of scattering, with a pulse-broadening time of 
0.8 ± 0.4 yits that is the lowest estimated yet towards the 
Crab from observations so far. Further, the pulsar's DM 
is determined to be 56.751 ± 0.001 pc cm~ 3 , which is sig- 
nificantly lower than those measured near the epochs of 
previous scattering measurements. 

Our measurements of DM and scattering, together with 
published data and the Jodrell Bank monthly cphcmeris, 
unveil a systematic and slow decrease in the Crab's DM 
and scattering over the past 10 yr. These variations are too 
large and smooth to be caused by the intervening ISM but 
can be attributed to the material within the nebula. Our 
analysis hints at there being large-scale inhomogeneities 
in the distribution of small-scale density structures in the 
nebular region, with a plausible interpretation involving 
many (~ 100) dense (~100 cm~ 3 ) structures. A variation 
in their number density or size can potentially lead to the 
observed changes in DM and scattering. Such a possibility 
can be further investigated by obtaining independent con- 
straints on the nebular electron densities (e.g. via free-free 
absorptions at low radio frequencies) and through future 
observations to monitor the pulsar's DM and scattering. 
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